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Using geometry to manipulate long-range correlation of light inside disordered media 
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We demonstrate experimentally that long-range intensity correlation for light propagating inside 
random photonic waveguides can be modified by changing the shape of the waveguide. The functional 
form of spatial correlation is no longer universal in the regime of diffusive transport and becomes 
shape-dependent due to the non-local nature of wave propagation. The spatial dependence of the 
correlation may be asymmetric for light incident from opposite ends of the waveguide. This work 
opens the door to control non-local effects in mesoscopic transport of waves by manipulating the 
geometry of random systems. 
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The diffusion model has been widely utilized to de¬ 
scribe wave propagation in disordered media, e.g., light in 
biological tissues, ultrasonic waves through cracked met¬ 
als, and electron wavefunctions in disordered conductors. 
It, however, ignores the interference of scattered waves, 
which lead to many prominent phenomena including An¬ 
derson localization, universal conductance fluctuations, 
and enhanced backscattering m- Extensive theoreti¬ 
cal and experimental studies in the past three decades 
have illustrated that mesoscopic transport of both classi¬ 
cal and quantum mechanical waves is governed by wave 
interference effects HIS]. 

One important consequence of wave interferences in 
random media is the correlations in the fluctuations of 
scattered intensities laizi- The interferences between 
waves scattered along independent paths give rise to in¬ 
tensity correlation on the scale of wavelength, one crossing 
of paths generates long-range correlation beyond the mean 
free path, and two crossings leads to an infinite-range cor¬ 
relation [SE]. The non-local correlations have a direct 
consequence for the coherent control of light transmission 
through random media via wavefront shaping cni, which 
has advanced rapidly in the past few years due to poten¬ 
tial applications to deep tissue imaging [IIHS]- Indeed 
focusing light to a single speckle simultaneously bright¬ 
ens nearby speckles, and hence reducing the contrast of 
focusing [I4j . However, the long-range correlation also 
facilitates the enhancement of total transmission m with 
an optimized wavefront despite of the limited degree of 
input control m- Moreover, the spatial correlation of 
intensity inside the random medium [TTUST)] affects focus¬ 
ing and energy deposition into the sample |2Ij . Therefore, 
manipulating the non-local correlation can open up a new 
avenue to controlling waves inside random media. 

Typically the magnitude of long-range correlation is 
small, but it becomes significant in strongly scatter¬ 
ing media, especially when localization regime 0 is ap¬ 
proached O El [221125] . Experimentally long-range corre¬ 
lation have been observed not only in space, but also in 
time, frequency, angle, and polarization, but most mea¬ 


surements are performed on transmitted or reflected light, 
i.e. outside the random media [T51I2SH55] . Modifications 
of the correlations of transmitted light have been realized 
with two techniques: (i) varying the probe beam spot 
size on a wide disordered slab HIEIETI, and (ii) in¬ 
serting a constriction, e.g., a pin hole, inside a random 
medium |291I34| . However, the possibility of manipulating 
long-range correlation inside the random media has not 
been explored. This is at least in part due to experimental 
challenge of gaining a noninvasive access to the interior 
of the medium where the light propagates. 

We recently fabricated quasi-two-dimensional random 
waveguides to probe the transport inside from the third di¬ 
mension EniEaED- This experimental setup has enabled 
us to monitor directly how the long-range spatial correla¬ 
tion builds up inside the diffusive system [20]. Moreover, 
by reducing (or increasing) the width of a rectangular 
waveguide, we were able to enhance (or suppress) the 
crossing probabilities of scattering paths throughout the 
system and, therefore, to modify the magnitude of long- 
range correlation. However, the functional form of corre¬ 
lation remained unchanged, as it is universal for diffusive 
quasi-one-dimensional waveguides HZHH]. 

In this Letter, we experimentally demonstrate an ef¬ 
fective approach to manipulating the spatial dependence 
of long-range intensity correlation inside random media. 
This is accomplished by fabricating photonic waveguides 
with cross-section varying along their length. The func¬ 
tional form of the long-range correlation is modified inside 
waveguides of different shapes because the crossing prob¬ 
ability of scattering paths is affected non-uniformly in 
space. Our approach enables global optimization via sys¬ 
tem geometry and it is applicable to other types of waves 
such as acoustic waves and matter waves. Besides the 
fundamental importance, manipulating the long-range 
correlation of waves inside random systems is useful for 
imaging and focusing into multiply scattering media us¬ 
ing wavefront shaping uniEiiisT] because it affects such 
aspects as focusing contrast, degree of control, as well 
as energy deposition inside the medium. Therefore, our 
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approach can provide an additional degree of freedom for 
controlling wave transport in scattering media. 

To illustrate the effects of waveguide geometry on long- 
range spatial correlation, we start with a theoretical anal¬ 
ysis of two-dimensional (2D) waveguides. The reflecting 
sidewalls confine the light inside the waveguide where 
scattering and diffusion take place within r = [y, z) plane 
with z being the axial direction. Light transport in the 
random waveguide is diffusive, and the non-local intensity 
correlation are dominated by the long-range correlation 
C 2 [Sld^]. The 2D correlation function C' 2 (ri; r 2 ) between 
two points ri = (?/i,zi) and r 2 = {y 2 ,Z 2 ) is calculated 
with the Langevin approach [13 HZl EHHini , see [H] for 
details. 

Let us consider the simplest case of linear tapering, 
namely, the waveguide width W(z) increases or decreases 
linearly along the waveguide axis z. Figure 1 shows the 
magnitude of C 2 , C 2 (r; r), in three waveguides with W{z) 
being constant (a), linear increasing (b) or linear decreas¬ 
ing (c). The 2D distributions across the waveguide are 
clearly different in the three cases, revealing the waveg¬ 
uide geometry has a significant impact on the growth of 
C 2 . In Fig. l(d-f), the correlation functions C 2 (z:i;z 2 ) of 
the cross-section averaged intensity [H] further illustrates 
the difference: in the waveguide of increasing W{z), the 
correlation function stay nearly constant for most values 
of zi and Z 2 , while in the waveguide of decreasing width, 
the correlation function exhibits more rapid variation over 
zi and Z 2 . These results suggest that the range of spatial 
correlation is increased (or decreased) in the gradually 
expanding (or contracting) waveguide, as compared to 
the waveguide of constant width. 

For a more quantitative comparison, the magnitude of 
C 2 of the cross-section averaged intensity, i.e., C 2 {z;z), 
is plotted in Fig. 2(a) for six waveguides of same length 
but different geometry. To compare the shape of these 
curves, the maximum value of each curve is set to 1. After 
the normalization, the two curves for the constant widths 
of 10 /im and 60 /rm coincide and agree to the universal 
functional form. In the expanding waveguide, C' 2 (z;z) 
increases more rapidly at the beginning and levels off 
when light diffuses deeper into the waveguide. This is 
attributed to the higher crossing probability of scattering 
paths near the front end of the waveguide where the cross 
section is narrower. As the width increases with z, the 
crossing probability is reduced, and C 2 rises more slowly. 
The contracting waveguide exhibits the opposite trend: 
the magnitude of C 2 grows more quickly in the second 
half of the waveguide due to enhanced crossing probability. 
The larger the tapering of the waveguide cross section, 
the bigger the change in the spatial dependence of € 2 - 

Figure 2(b) plots the correlation function C 2 {z-,L) for 
two points z and L of cross-section averaged intensity 
of the six waveguides studied above. After normalizing 
the maximum value to 1, C 2 {z; L) for the constant-width 
waveguides coincide, whereas in the expanding waveguide. 
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FIG. 1: (Color online) Calculated spatial long-range intensity 
correlation for the constant-width and two types of tapered 
2D random waveguides. The waveguide length L = 80 /rm, the 
transport mean free path I = 2.2 y,m, and the diffusive absorp¬ 
tion length = 26 pm. The waveguide in (a,d) has a constant 
width W = 10 pm; in (b,e) W{z) increases linearly from 10 
pm to 60 pm, while in (c,f) W{z) decreases linearly from 60 
pm to 10 pm. (a-c) show spatial distribution of the magnitude 
of long-range correlation function, C' 2 (r;r) for three geome¬ 
tries. (d-f) show long-range correlation function < 72 ( 21 ; 22 ) of 
cross-section averaged intensity m for the same geometries. 
The maximum value is normalized to 1 for comparison. The 
differences in these plots reveal that the waveguide geometry 
has a signihcant impact on the magnitude and range of (72. 
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FIG. 2: (Color online) Comparison of calculated long-range 
correlation in six waveguides with different degrees of taper: 
two with constant widths of 10 pm (solid black line) and 60 pm 
(dash-dotted magenta line); two with width linearly increasing 
from 10 pm (thick dashed blue line) or 20 pm (thin dashed 
blue line) to 60 pm; and two with width linearly decreasing 
from 60 pm to 10 pm (thick dotted red line) or 20 pm (thin 
dotted red line). Other parameters are the same as in Fig. 1. 
Both < 72 ( 2 ; 2 ) (a) and < 72 ( 2 ; L) (b) clearly demonstrate that 
while the functional form of long-range correlation is universal 
for uniform waveguides, it is strongly modified in the tapered 
ones. 


it increases more rapidly with z in the first half of the 
waveguide, and then increases at a much lower rate in the 
second half. In the waveguide with decreasing width, the 
behavior is opposite. 

These results confirm that the functional form of long- 
range correlation is modified by tapering of the waveguide. 
The deviation from the universal functional form in the 
constant-width waveguides is larger when the degree of 
tapering is stronger. The difference in the correlation 
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functions in expanding and contracting waveguides re¬ 
veals that (72 (ri; r 2 ) is no longer symmetric because one 
waveguide is a mirror image of the other. In other words, 
the long-range intensity correlation function for light in¬ 
put from one end of the tapered waveguide is different 
from that with input from the other end. This behavior is 
distinct from that of the constant-width waveguide whose 
two ends are equivalent. 

Next, we conduct the experiments. The 2D disordered 
waveguides are fabricated in a silicon-on-insulator (SOI) 
wafer with a 220 nm silicon layer on top of a 3 fj,m buried 
oxide. The patterns are written by electron beam lithog¬ 
raphy and etched in an inductively-coupled-plasma (ICP) 
reactive-ion-etcher (RIE). Each waveguide contain a 2D 
random array of air holes that scatter light. The air hole 
diameters are 100 nm and the average (center-to-center) 
distance of adjacent holes is 390 nm. The waveguide 
walls are made of triangle lattice of air holes (the lat¬ 
tice constant of 440 nm, the hole radius of 154 nm) that 
has a complete 2D photonic bandgap for the in-plane 
confinement of light. 

The monochromatic beam from a tunable CW laser 
source (HP 8168E) is coupled into the empty waveguide 
by an objective lens of numerical aperture (NA) 0.4. The 
light is transverse-electric (TE) polarized, i.e., the electric 
field is in the plane of the waveguide. After propagating 
through the empty waveguide, the light is incident onto 
the random array of air holes inside the waveguide. The 
front end of the random array is uniformly illuminated 
along the y direction. The light undergoes multiple scat¬ 
tering in the 2D plane of waveguide. Some of the light 
is scattered out of plane and imaged by a 50x objective 
lens (NA = 0.42) onto an InGaAs camera (Xeva 1.7-320). 

From the optical image, the spatial distribution of light 
intensity inside the waveguide I{y,z) is extracted. To 
smooth out the short-range fluctuations, I{y,z) is aver¬ 
aged over the cross-section of the waveguide to obtain 
the cross-section-averaged intensity Iv{z). The spatial in¬ 
tensity correlation C{zi,Z 2 ) is then computed from Iv{z). 
After removing the short-range contributions, C{zi, Z 2 ) is 
dominated by long-range correlation C 2 - The contribu¬ 
tion of (73, which is on the order of 1 / 5 ^, is negligible as 
g 1 in our waveguides. 

The relevant parameters for light transport in the dis¬ 
ordered waveguide are the transport mean free path I and 
the diffusive dissipation length ■ The dissipation results 
from out-of-plane scattering, which can be treated simi¬ 
larly as absorption |35j . From the disordered waveguides 
with constant width, we find = 26 ym and £ = 2.2 ym 
by fitting the measured I{z) and C{zi,Z 2 = zi) [TT]. The 
tapered waveguides have the same density and diameter 
of air holes, thus the values of and £ ym are identical. 

Figure 3(a,b) are the scanning electron microscope 
(SEM) images of an expanding waveguide and a contract¬ 
ing waveguide. The measured correlation functions for 
the cross-section averaged intensity inside the two waveg¬ 


uides, C{zi = z,Z 2 = L), are shown in Figure 3(c). To 
compare the functional form, the measured C{z,L) is 
normalized. The experimental data clearly show that the 
dependence of C{z, L) on z is very different for the two 
tapered waveguides, which agree well to the calculation 
results. 

Since the waveguide geometry in Fig. 3(b) is the mirror 
image of the one in Fig. 3(a), the (7(z, L) for light input 
from the left end of the former is equivalent to that 
with input from the right end of the latter. As (7 is 
dominated by long-range correlation function, this result 
implies C 2 becomes asymmetric. Therefore, by tapering 
the waveguides we can modify the crossing probabilities 
of scattering paths and vary the long-range correlation. 
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FIG. 3: Experimental measurement of long-range intensity 
correlation inside the tapered waveguides. (a,b) Top-view 
SEM images of fabricated quasi-2D disordered waveguides with 
linearly increasing (a) or decreasing (b) width. The width of 
waveguide in (a) increases from 10 ym to 60 ym, and in (b) it is 
opposite. Both have the same length L = 80 ym. (c) Magnified 
SEM images showing the air holes distributed randomly in the 
tapered section of the waveguide and the triangle lattice of 
air holes in the reflecting sidewalls, (d) Measured long-range 
correlation function for the cross-section-averaged intensity 
C{z,L) inside the tapered waveguides shown in (a) and (b). 
The blue circles (green squares) represent experimental data 
for the waveguides with increasing (decreasing) width, and 
the dashed lines are theoretical results. 


Finally, we modulate the waveguide cross section in a 
non-monotonic manner for further manipulation of long- 
range intensity correlation inside the random waveguide. 
The waveguide shown in Fig. 4(a) has the width W in¬ 
creasing linearly in the first half of the waveguide and 
decreasing in the second half. This geometry, unlike the 
tapered waveguides studied above, is symmetric with re¬ 
spect to the center (z = T/2), thus the spatial intensity 
correlation function is the same for light incident from 
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either end of the waveguide. Figure 4(b) shows the spa¬ 
tial distribution of light intensity inside the waveguide 
with input from the left end. The short-range intensity 
fluctuations seen in Fig. 4(b) are smoothed out after 
averaging the intensity over the cross section, leaving only 
the long-range contributions to the intensity correlation 
function C(zi,Z 2 ). Figure 4(c) plots C{z,L), which in¬ 
creases initially at a slow rate as z approaches T/2, but 
turns into a sharp rise once z passes L/2 and approaches 
L. This is because the crossing probability of scattering 
paths is first reduced as the waveguide is expanding in 
z < T/2, and then enhanced in z > L/2 as the cross 
section decreases. Therefore, the crossing probability can 
be controlled by modulating the waveguide width, which 
changes the spatial dependence of long-range correlation 
function. Figure 4(d) shows the intensity correlation func¬ 
tion C(zi = z, Z 2 = L/2). It first increases monotonically 
as z moves from 0 to L/2, and then decreases slightly for 
z from L/2 to L. The experimental data (solid circles) 
are in good agreement to the theoretical results (dashed 
lines) in Fig. 4(c,d). 
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FIG. 4: Long-range correlation in a quasi-2D disordered waveg¬ 
uide whose width varies non-monotically. (a) Top-view SEM 
image showing the waveguide width W increases linearly from 
10 /rm at z = 0 to 60 fim at z = 40 /rm and then reduces 
linearly down to 10 fim at z = 80 pm. Other structnral param¬ 
eters are the same as the waveguides in Fig. 3. (b) An optical 
image of the intensity of scattered light from the disordered 
waveguide. The wavelength of the probe light is 1510 nm. (c) 
Long-range correlation function C{z, L) for the cross-section 
averaged intensities at z and L in the waveguide shown in (a). 
C(z, L) displays a sharp change in the growth rate before and 
after z passes L/2. (d) Long-range correlation function C{z, L) 
for the cross-section averaged intensities at z and L/2 in the 
waveguide shown in (a). C{z,Ll2) increases monotonically in 
the first half of the waveguide and decreases slightly in the 
second half. In (c,d), solid circles represent experimental data 
and the dashed curves are obtained by numerical calculation. 


In summary, we have demonstrated an efficient control 
of long-range intensity correlation inside random media by 


tailoring the shape of a random system without modifying 
the structural disorder. By modulating the cross section 
of a disordered waveguide, we manipulate the crossings of 
the scattering paths inside the waveguide and change the 
functional form of long-range correlation from the univer¬ 
sal one for constant-width waveguides. In addition, the 
long-range correlation function can be made asymmetric 
with the spatial dependence different for light incident 
from one end or the other end of the waveguide. This 
work opens a new direction for coherent control of wave 
transport in random media. 
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